Response delay as a strategy for survival in fluctuating environment 
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Response time-delay is an ubiquitous phenomenon in biological systems. Here we use a simple 
stochastic population model with time-delayed switching-rate conversion to quantitatively study the 
biological influence of the response time-delay on the survival fitness of cells living in periodically 
fluctuating and stochastic environment, respectively. Our calculation and simulation show that 
for the cells having a slow rate transition into fit phenotype and a fast rate transition into unfit 
phenotype, the response time-delay can always enhance their fitness during the environment change. 
Particularly, in the periodic or stochastic environment with small variance, the optimal fitness 
achieved by these cells is superior to that of cells with reverse switching rates even if the latter 
exhibits rapid response. These results suggest that the response time delay may be utilized by cells 
to enhance their adaptation to the fluctuating environment. 

PACS numbers; 



I. INTRODUCTION 

Most of living species are exposed to random fluctu- 
ations in their environment in the Nature. Environ- 
ment change can either be periodic, such as day-night 
interchanges, or be stochastic, such as nutrition sup- 
plements in random times. Thus how can a life-form 
adapt to the fluctuating environment to optimize its fit- 
ness within is critical for the species to survive in the 
nature [U, Q • Recently, with the progresses in the study 
of the gene regulation networks, the stochastic nature of 
the gene expression is revealed [3,13,0 , which is soon con- 
sidered to further understand optimized fitness of cells. 
Stochastic gene expression brings the isogenic popula- 
tion to distinct phenotyp^s which is known as phcno- 
typic heterogeneity 0, [1]. Phenotypic heterogeneity is 
a useful strategy for species to adapt to the fluctuating 
environment With this concept, various phenomena, 
including the bacterial persistence d, [l0| and opti- 
mized growth rate for the cells living in the fluctuating 
environment [ll|, [H, [H, [3, [Hi, were studied. These 
authors also considered comparison between the strategy 
of the heterogeneity in which the switching rate from the 
fit state to the unfit states in an environment of specific 
parameters is nonzero and the homogeneity in which the 
switching rate mentioned above is zero [nj, [H, the ef- 
fects of the static (switching rates do not change with 
environment) and dynamic (switching rates change with 
environment) transition between states [l^], and the in- 
fluence of the switching rates on the environment chang- 
ing rate [l^ , which are all important for the understand- 
ing of the survival of species on earth. 

Besides the stochastic gene expression, another impor- 
tant feature of gene regulation networks is time-delay. 
A time delay usually results from regulated transcrip- 
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tion and translation, and leads to behaviors quite dif- 
ferent from those without [l6j . In more complex gene 
networks, some motifs, such as feed-forward loops and 
gene cascades, induce time-delays of length even longer 
than a cell cycle [13, [H, [l!]. These time-delays have 
important effects on gene regulatory networks, and there 
were some assumption that cells can benefit from their 
time-delay motif in stochastically fluctuating environ- 
ments [20[. However, so far as we know, there is yet no 
work on the influence of the response time-delay to the 
fitness of cells living in the fluctuating environment. In 
this paper, we construct a stochastic population model 
with time-delayed switching-rate conversation to quan- 
tify the biological influence of the response time-delay on 
the survival of cells in fluctuating environments. 



II. MODEL 

To study the fitness of cell population, we present a 
stochastic population model [ll|, in which environment 
parameters fiuctuate between two different states Ea and 
Eb, and accordingly in which cells display two different 
phenotypes a and b. Cells of phenotype a are fit for 
growth in the environment Ea with fast growth rate 7f, 
while those of b is unfit for growth in Ea with slow 7u. 
Situation reverses in the environment E^ as shown in fig- 
ure 1(a). The environment compensation for the growth 
of the fit cells is A7 = 7f — 7u > 0. In addition, in 
these two environment states, cells stochastically switch 
between the two distinct phenotypes with switching rate 
kf into fit phenotype and fcu into unfit one. Some authors 
generally treated the switching-rate conversion simulta- 
neous with the environment change [Til . [l5| . But here in 
order to consider the influence of response time-delay, we 
set the switching rates conversion with a delayed time r 
after an state change; see figure 1(a). The delayed time 
T is a consequence of the gene networks structure when 
responding to environment changes. In this model, the 
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dynamics of the cell numbers in the fit (rif) and unfit 
(tt-u) states can be described as 



—Uf = 'YfUf + k^n^ — kfUf, < t < t; 
at 

d 

-rnn = 7u"-u + hnt ~ kuU^, < t < t; 
at 

-^rif = 7irif + fcfWu - fcu^f , T <t <T; 



(a) 



-^riu = 7u"-u + Knf - fef^u, T < t <T, 



(1) 



where T is the time cells dwell in an environment state 
before change to the other state. We assume t < T, 
which means that dynamic transition condition holds 
during the switching-rates conversion. When t > T, 
switching rates do not convert in this environment state 
which in fact becomes the static transition. However, 
static transition can also be a possible choice for cells, 
and below we consider above situation together. 

To investigate the fitness of cells, we consider 
the time-averaged population growth-rate (7) = 
(7u/u(<) + 7f/f(i)) = 7u + A7(/f(0), where ^t) = 
n-i/ (riu + n±), i = u, f is the fraction of cells in the unfit 
and fit phenotypes in the two environment states respec- 
tively. We assume in this model that the fractions of the 
cells in fit and unfit phenotypes convert when the cells 
respond to environment changes, that's, if the fraction of 
the cells in the fit phenotype is ft just before the envi- 
ronment turnover, it will he 1 — f± after. Hereafter we 
focus on the time evolution of the /f and its time-average 
((/f)) which accounts for the fitness of the cell popula- 
tion. By applying equation [TJ the dynamics of /f can be 
described by the following equations: 



jjf{t) = + (A7 -k^- kt).fi - A7/2, < i < r; 

jjf{t) = h + (A7 -k^- kt)ft - A7/2, r < t < T, (2) 

In different fluctuating environment conditions, T and r 
can be tuned, we solve equations [51 then average /f with 
time to evaluate the influence of response time delay on 
the survival fitness of cells. 



III. RESULT AND DISCUSSION 

Figured] (b) shows two examples of fi{t) changing with 
time in periodic and stochastic environments. In the pe- 
riodic environment state, cells spends a specific time T in 
each states, while in the stochastic case, the time spent 
in each states has a mean value T. Without losing gen- 
erality, we set r = 1. From the figure we can see that 
after considering the time delay, in the periodic envi- 
ronment there will be an extra turning point at time r 
after each environment change, while in the stochastic 
environment, only the environment states with duration 
time longer than r show a state turning point. Those 
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FIG. 1: (a) The scheme of the response time delay model 
of the fluctuating environment, (b) /(f) vs time in the peri- 
odic and stochastic environments, respectively. Correspond- 
ing fraction of the fit cells (/t) is plotted, while parameters 
are set as fcn — 0.5, ki = 1.5, A7 = 1 and r = 0.5. (c) 
Growth in a symmetrically periodic environment. The frac- 
tion of the fit cells (ft) is plotted with time for a period. In 
the figure we show the results of three specific r, 0, 0.5 and 
1. The parameters in the upper and lower plots are set as 
fcu = 1.5, fcf = 0.5, A7 = 1 and feu = 0.5, fef = 1.5, A7 = 1, 
respectively. 



brief environment states with duration time shorter than 
T are filtered as a result of phenotype switching system 
of cells. The additional turning points of the /f (t) in dif- 
ferent environments and filter out of the short duration 
environment states in the stochastic environment may in- 
troduce intriguing behavior in population fitness of cells 
with response time-delay, as we show below. 

To investigate the infiuence of response time-delay, we 
first consider a simple symmetrically periodic environ- 
ment with a specific time-delay [2l|. For different r, fi{t) 
reaches its periodic state as shown in figure 1(c). Two 
cases are shown there with inter-conversed switching- 
rates. For r = 0, the model returns to what is studied 
by Thattai et. al, and so does when r = 1, in which k-a 
denotes the rate switching to the fit state and fcf to the 
unfit one. Therefore, from figure [1] (c), we can see that 
the pattern of ff{t) in the upper graph when r = is 
the same as the lower one for t ~ I, and so does when 
T = 1 and r = . Since (/f(t)) reflects the fitness, we 
then calculate its averaged value. For these two delayed 
time discussed above, one can easily notice that the case 
with larger switching rate to the fit phenotype generally 
have higher (ft); see figure 2(a). When < t < 1, the 
tuning point in ff{t) (figure [1] (b), (c)) makes (/f) be- 
haviors quite different from those without time delays. 
After some mathematical analysis and with the help of 
numerical solution, we find that (ft) both decrease where 
T — > 0"*" and 1~, indicating a maximum at an intcrme- 
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(a) 



(b) 





FIG. 2: (Color online). Influence of the response time-delay 
on the time-averaged fraction of cells in fit phenotype. (a) 
Time-averaged /f is plotted with r. The solid and dash lines 
denote the case when fcn = 1.5 > kt = 0.5 and fcf — 1.5 > 
fcu — 0.5, respectively, as A7 = 1. Two cases of specific time- 
delays are also shown in the figure, (b) Parameter influence 
on the time-averaged ft. Here we set fcn — 1.5, A7 = 1, and 
kf < fcn from to 1.5. Insert shows kt > fcu from 1.5 to 5. 



diate r when /cu > kf. While kf > ku, fitness increases 
both as r — > 0"*" and 1~, indicating a minimum at an 
intermediate r (figure 2). Besides, there is another case 
when r > 1, when switching rates have no enough time 
to response to the environment change in neither envi- 
ronment states, which is just the case of static transition. 
In this case, the switching rate to the fit phenotype is half 
kf and half fcu in a whole period, and the average effect 
should be some combination of the k^ and kf. So the 
value of (/f } is between the two cases of t = and r = 1 
(figure 2(a)). 

The results shown above correspond to certain sets of 
parameters. In all there are three parameters that can 
be tuned in this model. A7 represents the growth com- 
pensation in different environment. Large A7 indicates 
much faster cells growth in fit environment than in unfit 
one, thus it has significant effects on (/f), but limited 
effects on the trend of (/f) to the time-delay r. Larger 
A7 induces higher (/f), but make smaller the difference 
of (/f ) induced by different t. However, an extremum at 
certain r always exist (data not shown). The other two 
tunable parameters are k^ and kf. They show similar 
effects in our modcl(cquation 1). Here wc first keep k^ 
fixed and vary kf to see its influence on the (/f). As 
shown in figure 2(b), when kf = fcu, exchange of switch- 
ing rates become trivial in the model, thus delay times 
have no influence on (/f) and (/f) remains constant at 
different t. When kf < k^, there is a maximum (/f) at 
a certain r, even in the homogenetic case when {kf = 0). 
In addition, in the delayed time (0~ r with fcu), larger 
k-a/kf indicates a faster transition into the fit phenotype, 
thus results in a steeper slope of the (/f) while reaching 



its maximum, which means more significant influence of 
response time-delay on increasing the fitness in case of 
larger k^/kf (figure 2(b)). A contrary effect is observed 
for kf > fcu, which always shows a minimum(figure 2(b) 
insert). The slope of the decreasing (/f) is steeper with 
larger fcf/fcu, which means the response time-delays also 
show more significant negative influence on the fitness 
with larger kf/k^ (figure 2(b) insert). These results ver- 
ify that the influence of the response time delay on the 
fitness of cells is universal and does not strictly depend 
on state parameters. 

These results show some interesting behavior of the re- 
sponse time delay on the fitness of cells. First, if kf > k^ 
which means a larger transition rate into fit phenotype in 
each environment, the best strategy to enhance fitness is 
to reduce the response time delay, as indicated in figure [2] 
(b) and larger fcf/fcu introduces a more severe defect in 
fitness. Second, above results tells us that cells do not 
always need to adopt kf > fcu to achieve a higher fitness. 
If cells use some special strategy to reverse the value of 
the two switching rates to set k^ > kf, the response time 
delay can always approve cells' fitness. Besides, there is a 
specific delay time ti after which (/f ) is larger than what 
cells can reach without time delay in the similar set of 
parameters; see figure 2(a). In the region ti < t < 1, 
the response time delays prove to be an effective way for 
cells of optimizing their fitness. In addition, considering 
the constant presence of an intrinsic response time-delay 
in gene networks, there are two switching rates cells can 
choose. For example, when the delay time is longer than 
T2, a strategy which sets fcu to be the larger one gets a 
better fitness (figure 2(a)). The parameter-dependence 
of Ti and T2 are shown in figure 3(a). The value of ti is 
always smaller than 1, which confirms the independence 
of parameters of these effects, relatively. Figure 3(a) also 
shows a large parameter space to be realized to optimize 
fitness by response time-delay. Third, the static transi- 
tion strategy (t > 1) can be a choice for cells in some 
special cases. For example, when the systemic time de- 
lay is large than a certain value with kf > fcu or shorter 
with kf < fcu (figure 2(a)), static transition is better than 
dynamic one in enhancing survival fitness, while in other 
conditions with the same switching rates, dynamic tran- 
sition strategy is still an better option. Finally, following 
the finding that ff{t) for t = and 1 are the same while 
exchanging kf and fcu, wc find that after a similar ma- 
nipulation of T with 1 — T and kf with fcu in equation 
1, this new equation becomes one describing a system 
with a pre- response time t. Therefore, the (ff) depen- 
dence on the pre-response time r in the new equation for 
kf > fcu is in fact the same as those with time-delay r 
in the former equation with kf < fcu, which means that 
an advance-response may induce a better fitness for the 
population when the cells adopt a faster switching rate 
to the fit phenotype [2^ . 

Then we consider a simple stochastic environment 
whose duration-time distribution is described by T = 
'^i^i Ti, where Ti obeys an exponential distribution with 
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FIG. 3: (a)Parameter-dependence of the two special delay 
times. Black and grey lines denote ri and T2, respectively. On 
the left is their dependence on A7, with kf = 10, fcu = (solid 
line), kf = 5, k^ — 0.1 (dashed line), and kf = 1.5, fcu = 1 
(dotted line). On the right shows their dependence on kf, 
with fcn — 1.5 and A7 = 1 (solid line), 10 (dashed line), and 
30 (dotted line), respectively, (b) and (c) Cell growth in the 
fluctuating environment. Time-averaged ff is plotted with 
T, with n = 1 in (b) and n = 10 in (c). The black and 
grey lines denote the case when fcu — 1.5 > kf = 0.5 and 
kf — 1.5 > fcu = 0.5, respectively. A7 = 1. 



mean value l/n. So the mean value of T is 1 and its vari- 
ance 1/n. When n = I, T itself obeys an exponential 
distribution, and when n — > 00, it returns to a determin- 
istic periodic environment state. A general form of ft{t) 
is shovifn in figure 1(b). Our results show, when n = I and 
fcu > fcf , time delays increase the fitness without reaching 
a maximum (figure 3(b)). Similarly, after converting the 
two switching rates to set < kf, (ft) monotonically 
decreases asymptotically to a value with increasing time- 
delay. However, in contrast to the former results, (/f ) for 
fcu < kf is always larger than that while k^ > kf. This 
phenomenon shows that in this stochastic environment 
the response time-delay is always a defect for the fitness, 
and tuning k^ < kf is simply a better strategy [20l | . When 
in cases of larger n, results gradually return to the for- 



mer ones. Results in case of n = 10 are shown in figure 
3(c), which shows almost all the characteristics we dis- 
cuss above. These results show that the variance may 
have significant influence on the cells' fltness just as the 
averaged value of the fluctuating environment descrip- 
tion. Large variance implies an existence of very short- 
or very long-duration environment states. Short ones are 
filtered by the cells' phenotype switching system as we 
discussed above, and this effect is similar with the static 
transition. For the long ones, there are enough time for 
cells to alter their phenotype with the prescribed switch- 
ing rates according to the environment change. There- 
fore, these two effects erase maximum or minimum (ff) 
induced by time delay, and (ff) monotonously converge 
to the static value with increasing r. Similar results hold 
in non-symmetric periodic environment description. If 
the ratio of the duration time of the two distinct envi- 
ronment states becomes larger, (/f ) of the k^ < kf also 
become larger and those of reversed parameters become 
smaller accordingly. Finally they separate and don't in- 
tersect just as what happens in the n = I stochastic 
environment (data not shown). 

In this work, we use a simple stochastic population 
model with time-delayed switching-rates conversion to 
study the influence of the response time-delay on the fit- 
ness of cell population. We find that in some cases such as 
fcu > fcf , time delays can optimize the population fitness, 
while in some other cases, delay time don't prove to be so 
positive. We also find that in the stochastically fluctuat- 
ing environment states, variance may be another impor- 
tant factor determining whether time delay can be an ad- 
vantage for the population fitness or not. However, since 
we only considered the exponential distribution of du- 
ration time above, other type of stochastic environment 
should also be considered so as to give a better under- 
standing of population fitness in such environment [20| . 
Besides, this model only uses the average growth-rate 
to account for population fitness while neglecting other 
factors such as energy consumption, thus other environ- 
ment state factors can be included to give a more com- 
prehensive understanding of the issue, and experimental 
verification of mentioned results is also needed. 

This work was supported in part by Tsinghua Basic Re- 
search Foundation and by the National Science Foun- 
dation of China under Grant No. 10704045 and No. 
10547002. 



[1] Schaffer, W.M., Am. Nat. 108, 783 (1974). 

[2] Mroz, I., P§kalski, A. and Sznajd-Weron, K., Phys. Rev. 

Lett. 76, 3025 (1996). 
[3] Kaern, M., Elston, T.C., Blake, W.J. and CoUins, J. J., 

Nature Rev. Genet. 6, 451 (2005). 
[4] Samoilov, M.S., Price, G. and Arkin, A.P., Sci. STKE 

2006, rel7 (2006). 
[5] Kaufmann, B.B. and van Oudenaarden, A., Curr. Opin. 



Genet. Dev. 17, 107 (2007). 
[6] Lu, T., Shen, T., Bennett, M.R., Wolynes, P.G. and 

Hasty, J., Proc. Natl. Acad. Sci. USA 104, 18982 (2007). 
[7] Sasaki, A. and EUner S., Evolution 49, 337 (1995). 
[8] Balaban, N.Q., Merrin, J., Chait, R., Kowalik, L. and 

Leiler, S., Science 305, 1622 (2004). 
[9] Kussell, E., Kishony, R., Balaban, N.Q. and Leibler, S., 

Genetics 169, 1807 (2005). 



5 



[10] Bishop, A.L., Rab, F.A., Sumner, E.R. and Avery, S.V., 

Mol. Microbiol. 63, 507 (2007). 
[11] Thattai, M. and van Oudenaarden, A., Genetics 167, 523 

(2004). 

[12] Kussell, E. and Leibler, S., Science 309, 2075 (2005). 
[13] Kussell, E. and Leibler, S., Phys. Rev. Lett. 97, 068101 
(2006). 

[14] Acar, M., Mettetal, J.T., and van Oudenaarden, A., Na- 
ture Genet. 40, 471 (2007). 

[15] Pang, N.N. and Tzeng, W.J., Bull. Math. Biol. 70, 210 
(2008). 

[16] Bratsun, D., Volfson, D., Tsimring, L.S. and Hasty, J., 
Proc. Natl. Acad. Sci. USA 102, 14593 (2005). 

[17] Mangan, S. and Alon, U., Proc. Natl. Acad. Sci. USA 
100, 11980 (2003). 



[18] Rosenfeld, N. and Alon, U., J. Mol. Biol. 329, 645 (2003). 
[19] Mangan, S., Zaslaver, A. and Alon, U., J. Mol. Biol. 334, 
197 (2003). 

[20] Dekel, E., Mangan, S. and Alon, U., Phys. Biol. 2, 81 
(2005). 

[21] We have simulated the effect of distinct response delay- 
time for two environment states, while little influence on 
the trend of the (/f ) dependence on r is showed, except 
the case when the larger r dorminates when the system 
becomes of a static transition. So here we only consider 
a specific reponse delay time. 

[22] Baliga, A.S., Science 320, 1297 (2008). 

[23] Tagkopoulos, I., Liu, Y.C. and Tavazoie, S., Science 320, 
1313 (2008). 



